Positron attachment to the He doubly excited states 
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The projection method is used to demonstrate the existence of positron attachment to three 
doubly excited states of helium. The e+He(2s^ ^8*=), e+He(3s^ ^8*=), and the e+He(2s2p ^P°) states 
have binding energies of 0.447 eV, 0.256 eV and 0.486 eV respectively. These energies were computed 
with the stochastic variational method and the conhguration interaction method. These states will 
exist as resonances in the e^-He continuum and the e^He(2s^ ^S"^) state could be detectable in the 
e"'"-|-He coUision spectrum. The existence of a series of e+He(ns'^ ^S^) resonances associated with 
the E[e(ns^) double Rydberg series is also predicted and an explicit calculation demonstrating the 
existence of the e"'"He(3s^ ^S*^) state is reported. 
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In this paper the ability of a positron to attach itself to 
the doubly excited states of helium is demonstrated by 
explicit calculation using the same projection operator 
approach that was originally used in some of the earliest 
calculations of the helium doubly excited spectrum [l|, ■ 
Besides the intrinsic interest in such exotic Coulomb sys- 
tems, the result provides a pathway to providing exper- 
imental confirmation that positrons can be attached to 
electrically neutral atoms to form bound states. 

It is now widely accepted that positrons can form 
bound states with a variety of atoms [^-Q- While the 
evidence for positron binding is very strong, it is mainly 
derived from calculation. Binding energies range from 
0.0129 eV in the case of e+Na Q to 0.50 eV for the 
e+Ca ground state 

There is solid experimental evidence that positrons can 
form bound states with a variety of molecules. The en- 
ergy resolved positron annihilation cross sections for a 
number of molecules (e.g. CaHg, C6II14) show features 
that have been identified as Feshbach resonances formed 
by the trapping of positrons to the vibrationally excited 
states [1] . Positron binding to vibrationally excited states 
of molecules is thought to be the mechanism responsi- 
ble for the large positron annihilation rates observed for 
many molecules in gas phase positron annihilation spec- 
troscopy experiments 

While the experimental evidence of positron binding 
to molecules is good, there is no experimental evidence 
that could be construed as demonstrating the existence 
of positron-atom bound states. One possible signature 
would be the existence of resonant structures associ- 
ated with atomic excited states in the positron scattering 
spectrum. Years of experimentation, however, have re- 
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vealed little evidence for the existence of resonant states 
in positron-atom scattering spectra [3. [lOl. [T]| . 

A number of schemes have been put forward to 
demonstrate the existence of positron-atom bound states 
[l^ - [T5| . The most recent proposal suggested that 
positron scattering experiments be performed on open 
shell transition-metal atoms having polarizabilities and 
ionization energies conducive to binding positrons [15|. 
Open shell systems are recommended since such systems 
would have low-lying excited states that could also bind 
a positron. Positron binding to these low-lying excited 
states would result in Feshbach resonances appearing in 
the low-energy annihilation cross section. However, the 
transition-metal atoms most likely to bind a positron rep- 
resent difficult propositions for experimentation. 

The present paper demonstrates that three of the 
doubly excited states of helium, namely the He(2s^ 
^S^), He(3s2 ^S^) and He(2s2p states can attach 
a positron with attachment energies exceeding 0.250 eV. 
The e+He(2s2 ^S^) and e+He(3s2 ^S^^) states should man- 
ifest themselves as broad resonances in the e^-l-IIe con- 
tinuum. A positron cannot excite the He(2s2p '^P°) state 
from the He(ls^ ^S'^) ground state since there is no ex- 
change interaction between the positron and electrons. 
These states can be regarded as the analogues of the 
triply excited negative ion resonances which are present 
in the electron-helium spectrum at 57 - 61 eV incident 
energy 

One motivation for the present investigation was the 
realization that the doubly excited states of helium have 
energetics very similar to those of the Mg atom. This 
binds a positron with a binding energy of 0.465 eV |29) . 
and also supports a prominent p-wave shape resonance 
in the elastic scattering channel at 0.096 eV incident en- 
ergy [29„ ^] . Table |I] compares energies of the doubly 
excited states of helium with the Mg(3£, 3£') states. The 
binding energy of the Mg+(3s) ground state is —0.55254 
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TABLE I: Comparisons of the energies of some He and Mg states. All energies are given in a.u. and given with respect to the 
He^+ and Mg^"'' thresholds. Three sets of helium energies are given. One set, Ecu are taken from complex rotation calculations, 
the two other sets are taken from projection operator calculations. The projection operator energies in the -Eqhq column come 
from calculations that use a Hylleraas basis, while those in the Eci column come from CI calculations as described in the text. 
There is no complex rotation energy for the He(2p^ '^P'^) state since it is a bound state and there is no state in the Mg spectrum 
that has been identified as being dominated by a 3p^ ^D*^ configuration. 



State 


-EcR 


-Eqhq 


Eci 


State 


-Ecxp [21] 


He+(2s) 


-0.500000 


-0.500000 


-0.500000 


Mg+(3s) 


-0.552536 


He(2s2 ^S") 


-0.777818 [22] 


-0.778774 [IJ 


-0.778781 


Mg(3s^) 


-0.833530 


He(2s2p ^P°) 


-0.760498 [23] 


-0.761492 [2] 


-0.761492 


Mg(3s3p ^P°) 


-0.733788 


He(2p2 ^P*=) 




-0.710500 [24] 


-0.710500 


Mg(3/ ^P") 


-0.569929 


He(2p2 ^D^) 


-0.701946 [25] 


-0.702817 [26] 


-0.702819 


Mg(3/ ^D^) 




He(2s2p ^P°) 


-0.69314 [22 


-0.692895 [28] 


-0.692897 


Mg(3s3p ^P°) 


-0.673824 


He(2p2 is-^) 


-0.621928 [22] 


-0.622744 [2] 


-0.622736 


Mg(3p2 iS=) 


-0.522446 



a.u. [21[ while the He"'' (2s) binding energy is —0.50 a.u. 
The binding energy of the Mg(3s^) ground state with re- 
spect to the Mg"''(3s) threshold is —0.2810 a.u, while the 
binding energy of the He(2s^) resonance with respect to 
the He+(2s) state is -0.2778 a.u. [12]. The respective 
dipole polarizabilities are 76.2 for He(2s^) and 71.3 Qq 
for Mg(3s2) 

The projection method ^, ^ provides a very simple 
calculation strategy for the identification of resonances 
and determination of their positions. In this method, 
the electrons are not allowed to occupy those low-lying 
states that could result in the auto-ionization of the sys- 
tem. The projection method energies, i^QHQ, of the he- 
lium doubly excited states in Table HI computed using 
Hylleraas basis sets, differ from those determined by the 
more dynamically complete complex rotation method by 
less than 0.001 a.u. The projection method has also been 
applied to the calculation of the positions of the He~ res- 
onances associated with the He doubly excited states [12] ■ 

The He"'' (Is) ground state is excluded from the con- 
figuration interaction (CI) wave function by Schmidt- 
orthogonalising the £ = single-particle electron orbital 
basis to the He+(ls) state. The Hamiltonian for the 
N — 2 electron and one positron system. 
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FIG. 1: Energy level diagram showing the positions of the 
He doubly excited states and the states with an attached 
positron. The position of the Ps(ls)-|-He^(2s) threshold is 
also shown. The axis on the right gives the positron collision 
energy (in eV) needed to excite these states. 
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(1) 



is then diagonalized in a large single particle basis of 
Laguerre type orbitals (LTOs). 

The CI method was first applied to the calculation of 
the He doubly excited states. The basis included 49 LTOs 
for £ = 0, and 50 LTOs for the other ^'s. The largest 
£ value used in these calculations was ^ = 8. The CI 
energies are given in Table HI and were extrapolated to 
the £ = oo limit using a procedure described shortly. 
They agree with the -Eqhq energies to within 10~^ a.u.. 



The e"'"He CI basis was constructed by letting the two 
electrons and the positron form all the possible configu- 
rations with the correct spin, parity and angular momen- 
tum, subject to the selection rules, 



max(£o,^i,^2) < J , 
min(£i,£2) < -^int , 

(_l)(^0+^l+^2) ^ _i ^ 



(2) 
(3) 
(4) 



In these rules £o, £i and £2 are respectively the orbital 
angular momenta of the positron and the two electrons. 
The number of LTOs for each £ was 15 with the exception 
of £ = 0, 1, 2, and 3 where 18, 18, 17, and 16 LTOs were 
used. The parameter Lint was set to 4. The largest £ 
in the orbital space was J = 12 for the e"'"He(2s^ '^S'^) 
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TABLE II: Calculated energies of some e'''IIe states. The 
CI calculations are also given with a J — > cxd correction as 
discussed in the text. The binding energies are denoted by e. 



State 




Metnoa 


T 

J 


{E)j (a.u.) 


e (a.u.) 


e (eV) 


e+He(2s2 ^ 


S") 


CI 


12 


-0.793537 


0.014756 


0.4015 






CI 


CXD 


-0.795058 


0.016277 


0.4429 






SVM 




-0.795190 


0.016409 


0.4465 


e+He(2s2p 




CI 


9 


-0.776306 


0.014814 


0.4031 






CI 


CXD 


-0.779362 


0.017869 


0.4863 


e+He(3s^ ^ 




CI 


12 


-0.468860 










CI 


CXD 


-0.481643 


0.009420 


0.2563 



state and J = 9 for the e+He(2s2p ^p") state. We define 
{E) J to be the energy of the calculation with a maximum 
single-particle orbital angular momentum of J . 

The main technical problem afflicting CI calculations 
of positron-atom interactions is the slow convergence of 
the energy with J 0, |3^. One way to determine the 
J — >■ oo energy, {E)oo, is to use an asymptotic analysis. 
It has been shown that successive increments, l^Ej = 
{E)j — {E)j-i, to the energy can written as an inverse 
power series |33l - [36| , viz. 



AEj 



Cp 



(J 



{J- 



(5) 



The J — > oo limit, has been determined by fitting sets 
of {E)j values to Eq. ([5|). The coefficients, Ae, Be and 
Ce are determined at a particular J from 4 successive 
energies {{E)j_3, {E)j_2, {E)j-i and {E)j). Once the 
coefficients have been determined it is easy to obtain the 
J ^ oo limit. Application of asymptotic series analysis 
to helium has resulted in CI calculations reprodu cing the 
ground state energy to an accuracy of 10^^ a.u. [36| . 

The CI energy of the e+He(2s2) state (see Table |TT| 
was —0.795058 a.u.. Subtracting this from the IIe(2s^) 
Eci of -0.778781 a.u. gives a binding energy of 0.016277 
a.u. This binding energy is an underestimate since the 
energy of the He(2s^) state in the CI basis used for the 
e+He calculation was —0.778771 a.u. The J ^ oo ex- 
trapolation contributed 10% to the binding energy. The 
e"''He(2s^) binding energy is 4% smaller than the bind- 
ing energy of the positron to the Mg(3s^) ground state, 
namely 0.01704 a.u. [2i|. 

The He(2s2p "^P") state also binds a positron with a 
binding energy of 0.017870 a.u. The surprisingly large 
binding energy is caused by the relatively small excita- 
tion energy of 0.051 a.u. from the IIe('^P°) state to the 
He(3p'=). This leads to the Re{^P°) state having a static 
dipole polarizability of 157.2 Aq. The e"'"IIe('^P°) positron 
binding energy is much larger than that of the e"'"Be('^P°) 
state which is only 0.000087 a.u ^S?]. 

The stochastic variational method (SVM) 0, [sl] was 
also used to determine the energy of the resonance state 
using the projection ansatz. The SVM diagonalizes the 
Hamiltonian in a basis of explicitly correlated gaussians 



(ECGs). The non-linear parameters of the ECG basis 
are optimized by a trial and error process. Such a pro- 
cess is possible since the ECG matrix elements of the 
Hamiltonian are very easy to compute. An orthogonal- 
izing pseudo-projector (OPP) [Sj, l39|, i4Q|] , is added to the 
Hamiltonian to exclude the He"*" (Is) state from being oc- 
cupied. The OPP operator is 



AP = A|ls)(ls| 



(6) 



Choosing A as a large positive number results in a ground 
state wave function which has a very small overlap with 
the He"*" (Is) state. The parameter A was set to 10® a.u. 
for the present calculations. The Is state was expanded 
as a linear combination of 12 gaussians. The SVM energy 
of the Hc(2s2) state was -0.778786 a.u., i.e. 5 x 10~® 
a.u. below the CI energy. The dimension of the largest 
SVM calculation of the e+He(2s2) state was 900 ECGs. 
The SVM binding energy of the positron to the He(2s^) 
state given in Table HIl was 0.016409 a.u. Examination of 
the convergence pattern suggests that the SVM energy is 
within 2 x lO""* a.u. of the variational limit. The SVM 
and CI binding energies for this state are in excellent 
agreement when the respective uncertainties arising from 
finite size basis sets are taken into consideration. 

The e"''He(2s^ ^S*^) system is also likely to support a 
^P° shape resonance just above the He(2s^ ^S'^) thresh- 
old. This is based on the similarity of the He and Mg 
polarizabilities and the positron attachment energies in 
the ^S*^ channel. The e^-Mg ^P° shape resonance was lo- 
cated at 0.00351 a.u. above the elastic scattering thresh- 
old and had a width of 0.00396 a.u. 

It is likely that there will be an infinite series of reso- 
nances associated with the set of He(ns^) doubly excited 
states. An investigation of the (m^"'", 2e^, e+) system 
revealed that this system remains bound when the mass 
TO^"*" — )■ [111 . Decreasing the m^"'" mass weakens the ef- 
fective strength of the m^+-e^ interaction and provides 
an analogue of the He^+-ns(e^) interaction. A first test 
was performed by a CI investigation of the e+Hc(3^,3^') 
systems. In this case the single particle basis was orthog- 
onalized to the He"*" (Is, 2s, 2p) states. The CI energy of 
the He(3s^) state is —0.354562 a.u. Since the removal en- 
ergy of the electron with respect to the He"*" (3£) thresh- 
old, —0.132340 a.u., is less than the positronium ground 
state energy of —0.25 a.u., the threshold for attaching a 
positron to the He(3s^) state is at —0.472222 a.u.. 

The CI calculation for the e"'"He(3s^) state gave an 
energy of —0.481643 a.u. The binding energy of this state 
is 0.009420 a.u. The stability of this system provides 
strong evidence for an infinite number of e+He(ns^) type 
resonances. It is likely that the rich resonance structures 
of the PsH system [42] will be replicated for positron 
interactions with the doubly excited helium atoms. 

Reference can be made to e~-|-He scattering experi- 
ments to give a first order estimate on the viability of 
experimental detection. A number of electron scattering 
experiments have demonstrated electron attachment to 
the He doubly excited states p^ - [20| . Experiments that 
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detect total cross sections involving ground state atoms 
and ions probably do not have a sufficiently large signal 
to background ratio to detect the e+He resonances. For 
example, He"*" ions were detected in the experiment of 
Quemener et al. [13] ■ There, the cross section for the 
creation of He+ varied by only 1% over the width of the 
He~(2s^2p) resonance. 

Higher signal to background ratios have been achieved 
in e~-He experiments that measure differential cross sec- 
tions. One example is an experiment which measured the 
energy of the ejected electron following auto-ionization 
[isj . Another is the measurement of the optical excita- 
tion functions associated with the formation of a number 
of He" excited states [l^ . 

The widths of the resonances and energy resolution of 
current positron beams is also an important considera- 
tion. Modern trap based positron beams can achieve an 
energy resolution of about 20-50 meV [1, H^. An ini- 
tial estimate of the widths can be made by equating the 



widths to those of their electron doubly excited parent 
states. The width of the He(2s2) state is F = 0.123 eV 
[13, the He(3s2) state is F = 0.082 eV [H, while the 



He(2s2p ^P°) state is F = 0.0081 eV [271 . The e+He(ns^) 



resonances should be wide enough to detect with current 
positron beam technology. The theoretically well-known 
positron-atom resonances are too narrow as in the case 
of hydrogen and sodium (i^ . l45j , involve atoms which do 
not naturally exist in gaseous form [2^ IH, lisl. o r require 
scanning over an energy range near 0.10 eV [23] which is 
difficult to achieve with existing positron beams. 
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